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Abstract. In this paper, we deduce a macroscopic strain gradient theory for plasticity from 
a model of discrete dislocations. 

We restrict our analysis to the case of a cylindrical symmetry for the crystal in exam, so 
that the mathematical formulation will involve a two dimensional variational problem. 

The dislocations are introduced as point topological defects of the strain fields, for which 
we compute the elastic energy stored outside the so called core region. We show that the 
T-limit as the core radius tends to zero and the number of dislocations tends to infinity of 
this energy (suitably rescaled), takes the form 

E= I (W(f3 e ) + v(Cur\ f3 e ))dx, 



Jci 

where /3 e represents the elastic part of the macroscopic strain, and Curl/3 e represents the 
geometrically necessary dislocation density. The plastic energy density tp is defined explicitly 
through an asymptotic cell formula, depending only on the elastic tensor and the class of 
the admissible Burgers vectors, accounting for the crystalline structure. It turns out to be 
positively 1-homogeneous, so that concentration on lines is permitted, accounting for the 
presence of pattern formations observed in crystals such as dislocation walls. 

Keywords : variational models, energy minimization, relaxation, plasticity, strain gradient 
theories, stress concentration, dislocations. 
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1. Introduction 

In the last decade the study of crystal defects such as dislocations has been increasingly active. 
The presence of dislocations (and their motion) is indeed considered the main mechanism of plastic 
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deformations in metals. Various phenomenological models have been proposed to account for plastic 
effects due to dislocations, such as the so called strain gradient theories. 

The aim of this paper is to provide a rigorous derivation of a strain gradient theory for plasticity 
as a mesoscopic limit of systems of discrete dislocations, which are introduced as point defects of the 
strain fields, for which we compute the elastic energy stored outside the so called core region. 

We focus on the stored energy, disregarding dissipation, and we restrict our analysis to the case of 
a cylindrical symmetry for the crystal in exam, so that the mathematical formulation will involve a 
two dimensional variational problem. 

Many theories of plasticity are framed within linearized elasticity. In classical linear elasticity a 
displacement of fl is a regular vector field u : fl — > R 2 . The equilibrium equations have the form 
DivC[e(u)] = 0, with C a linear operator from M 2x2 into itself, and e(u) := |(Vu + (Vu) T ) the 
infinitesimal strain tensor. The corresponding elastic energy is 

(1) f W{Vu)dx, 

Jn 

where 

(2) W{Cl = ±C£ : & 

is the elastic energy density, and the elasticity tensor C satisfies 

(3) ci|f ym | 2 < C£ : i < c 2 \e ym \ 2 for any ^eM 2 * 2 , 

where £ s y m := i(£ + £ T ) and c\ and C2 are two given positive constants. In this linear framework the 
presence of plastic deformations is classically modeled by the additive decomposition of the gradient 
of the displacement in the elastic strain f3 e and the plastic strain j3 p , i.e., Vu = f3 e + (3 P . The elastic 
energy induced by a given plastic strain is then 

(4) f W{\7u- (3 p )dx. 

Jn 

In view of their microscopic nature the presence of dislocations, responsible for the plastic deformation, 
can be modeled in the continuum by assigning the Curl of the field (3 P . The quantity Curl/3 P = fi is 
then called the Nye's dislocation density tensor, (see |19)V Inspired by this idea, the strain gradient 
theory for plasticity, first proposed by Fleck and Hucthinson [TU] and then developed by Gurtin 
(see for instance [13], P3]), assigns an additional phenomenological energy to the plastic deformation 
depending on the dislocation density, so that the stored energy looks like 

(5) f W{Vu - f3 p ) dx + [ p(Cur\ (3 P ) dx . 

Jn Jn 

Note that for a given preexisting strain (3 P , the minimum of the energy in ^ (possibly subject to 
external loads) depends only on Curl/3 P . Hence the relevant variable in this problem is a strain field 
(3 whose Curl is prescribed. The main issue with this model is the choice of the function ip in {5}. In 
fact the usual choice (see for instance jTS]) is to take ip quadratic, and then to fit the free parameters 
through experiments as well as simulations, as in [18] . This choice has the well-known disadvantage of 
preventing concentration of the dislocation density, that is instead allowed, for instance, by the choice 
<p(Curl/3) := |Curl/?|, recently proposed by Conti and Ortiz in [7] (see also [20]). 

The aim of this paper is to derive the strain gradient model ([5]) for the stored energy, starting from 
a basic model of discrete dislocations that accounts for the crystalline structure. As a consequence we 
will also get a formula for the function ip, determined only by the elasticity tensor C and the Burgers 
vectors of the crystal. 

A purely discrete description of a crystal in the presence of dislocations can be given after introduc- 
ing the discrete equivalent of (3 P and fi := Curl/3 P , following the approach of Ariza and Ortiz pQ. That 
model is the starting point for a microscopic description of dislocations. In the passage from discrete 
to continuum one can consider a intermediate (so to say) semi-discrete model known in literature as 
discrete dislocation model (see also [S], [5]), in which the atomic scale is introduced as an internal 
small parameter e referred to as the core radius. The gap between the purely discrete model and the 
discrete dislocation model can be usually filled through an interpolation procedure (see for instance 
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|21j). In this paper, to simplify matter, we will adopt the discrete dislocation model. 

In the discrete dislocation model a straight dislocation orthogonal to the cross section f2 is identified 
with a point io £ fi or, more precisely, with a small region surrounding the dislocation referred to 
as core region, i.e., a ball B s (xo), being the core radius e proportional to (the ratio between the 
dimensions of the crystal and) the underlying lattice spacing. The presence of the dislocation can be 
detected looking at the topological singularities of a continuum strain field (3, i.e., 

/ 13 -tds = b, 

JdB E (x ) 

where t is the tangent to dB £ (xo) and b is the Burgers vector of the given dislocation. 

We will focus our analysis to the case when there may be many dislocations crossing the domain fi. A 
generic distribution of N dislocations will be then identified with N points {xi}i t ... : N ( or equivalently 
with the corresponding core regions), each corresponding to some Burgers vectors bi belonging to 
a finite set of admissible Burgers vectors S, depending on the crystalline structure, e.g., for square 
crystals, up to a renormalization factor, S = {ei, e%, — e±, — e-z\. 

In our analysis we introduce a second small scale p £ 3> e (the hard core radius) at which a cluster 
of dislocations will be identified with a multiple dislocation. In mathematical terms this corresponds 
to introducing the span § of S on Z (i.e., the set of finite combinations of Burgers vectors with integer 
coefficients) and to representing a generic distribution of dislocations as a measure /i of the type 

N 
i=l 

where the distance between the xt's is at least 2p £ . The admissible strain fields (3 corresponding to p 
are defined outside the core region, namely in f2 \ UiB £ (xi), and satisfy 

(6) / P-td8 = £ i} for all i=l,...,N. 

It is well-known (and it will be discussed in detail in the sequel) that the plastic distortion due to the 
presence of dislocations decays as the inverse of the distance from the dislocations. This justifies the 
use of the linearized elastic energy outside the core region L>iB e (xi). On the other hand, considerations 
at a discrete level show that the elastic energy stored in the core region can be neglected. Therefore, 
the elastic energy corresponding to p and (3 is given by 

(7) E e (ji,0) := f W{P)dx, 

Jn e (fj,) 

where £l £ (p) = O \ UiB E (xi). By minimizing the elastic energy (O among all strains satisfying ([6]), we 
obtain the energy induced by the dislocation p. Note that, in view of the compatibility condition ([5]), 
this residual energy is positive whenever p, ^ 0. 

In the case W{(3) = |/3| 2 , with (3 a curl free vector field in f2\U,B E (^), this model coincides (setting 
(3 = V0, being the phase function) with the energy proposed by Bethuel, Brezis and Helein (see 
[4]) to deduce the so-called renormalized energy in the study of vortices in superconductors. In the 
context of dislocations this choice of W and (3 corresponds to a model for screw dislocations in an 
anti-planar setting. A T-convergence result in this framework has been studied in [21] in the energy 
regime E e fts | loge| that corresponds to a finite number of dislocations. 

In the planar context of edge dislocations, that we consider here, the connection with the Ginzburg- 
Landau model for vortices is more formal and the analysis has specific difficulties due to the vectorial 
nature of the circulation constraint © and the properties of the energy density W . A first analysis of 
this model has been perfomed by Cermelli and Leoni in [5] where they obtain an asymptotic expansion 
of the energy functionals E £ in l[7|). in the case of a fixed configuration of dislocations. In particular 
they deduce the corresponding renormalized energy representing the effect of the Peach-Koehler force. 

The purpose of this paper is to consider the energy functional E £ in jJJ , where the unknowns are 
the distribution of dislocations and the corresponding admissible strains, and to study its asymptotic 
behavior in terms of T-convergence as the lattice spacing e tends to zero and the number of dislocations 
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N £ goes to infinity. Under a suitable condition on the hard core scale, we can show that in the 
asymptotics the energy E £ can be decomposed into two effects: the self-energy, concentrated in the 
hard core region, and the interaction energy, diffused in the remaining part of Q. The general idea 
is that the T-limit of the (rescaled) energy functionals E e as e goes to zero is given by the sum of 
these two effects. For different regimes for the number of dislocations N e the corresponding T-limit 
will exhibit the dominance of one of the two effects: the self-energy is predominant for N £ <5< | loge|, 
while the interaction energy is predominant for N £ ^>> | loge|. In a critical regime N e « | loge| the two 
effects will be balanced and the structure of the limiting energy is the following 

(s, Jw^l^tH 

The first term is the elastic energy of the limiting rescaled strain and comes from the interaction 
energy. The second term represents the plastic energy and comes from the self-energy; it depends only 
on the rescaled dislocation density fi — Curl/3 through a positively homogeneous of degree 1 density 
function ip, defined by a suitable asymptotic cell problem formula. The constraint Curl (3 = /i comes 
from the admissibility condition (J6j) that the admissible strains have to satisfy at a discrete level. 

If (/3, /i) is a configuration that makes the continuous energy in (|8|) finite, we have necessarily that 
/3 sym belongs to L 2 (0; M 2x2 ), and has finite mass. A natural question arises: does \x belong to the 
space i? _1 (fi; M 2 )? We give a positive answer to this problem proving a Korn type inequality for 
fields whose curl is a prescribed measure, based on a fine estimate for elliptic systems with I^-data 
recently proved by Bourgain, Brezis and Van Schaftingen (see [3J, [5]). In particular we deduce that 
concentration of dislocations on lines is permitted, accounting for the presence of pattern formations 
observed in crystals such as dislocation walls, while concentration on points is not permitted. An 
additional feature of the limit energy is the anisotropy of the self-energy density inherited from the 
anisotropic elastic tensor and the class of the admissible Burgers vectors accounting for the crystalline 
structure. 

The idea of deducing continuous models by homogenization of lower dimensional singularities has 
been used in other contexts. Our result is very much in the spirit of earlier results on the asymptotic 
analysis for the Ginzburg-Landau model for superconductivity as the number of vortices goes to 
infinity, done by Jerrard and Soner in [16], and by Sandier and Serfaty in [23] (see also [24] and the 
references therein). 

As for the study of dislocations, a similar analysis was done by Garroni and Focardi Qllj) starting 
from a phase field model introduced by Koslowski, Cuitino and Ortiz ([17]) and inspired by the 
Peierls-Nabarro model (see also |12j). 

The plan of the paper is the following. In Section[2]we compute hcuristically the asymptotic behavior 
of the self and the interaction energy in terms of the number of dislocations N £ and the core radius 
e. In Section [3] we introduce rigorously the mathematical setting of the problem, defining the class 
of admissible configurations of dislocations, the class of the corresponding admissible strains and the 
rescaled energy functionals. In Section [4] we provide the asymptotic cell problem formulas defining the 
density of the plastic energy if. In Section [5] we prove a Korn type inequality for fields with prescribed 
curl . In Section [6] we give our main result concerning the T-convergence of the energy functionals in 
the critical regime N £ « | loge|, while Section [7] and Section [8] are devoted to the sub critical and to 
the super critical case respectively. 

2. Heuristic for the scaling 

In this section we identify the energy regimes, as the internal scale e tends to zero, of the elastic 
energy induced by a distribution of dislocations fi £ :— Yli=i^xi, where we omit the dependence 
on e of the x^s. The idea is that if the scaling of the energy of a given sequence of distributions of 
dislocations is assigned, this provides a bound for N e . In this section, in order to identify the relevant 
energy regimes, we do the converse: we assume to know the number of dislocations N £ present in the 
crystal, and we compute (heuristically) the corresponding energy behavior as e goes to zero. 

As already mentioned in the introduction, the general idea is that the energy will be always given by 
the sum of two terms: the self-energy, concentrated in a small region (the hard core region) surrounding 
the dislocations, and the interaction energy, diffused in the remaining part of the domain O. We will 
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first introduce and compute heuristically the asymptotic behavior of these two terms separately. Then 
we will identify the leading term (which depends on the behavior of N £ ) between these two energies. 

2.1. Self-energy. The self-energy of a distribution of dislocations fj, e := Yli^i £i$xi 1S essentially given 
by the sum of the energies that would be induced by a single dislocation /i* := (ii6 Xi , i = 1,...N S 
(i.e., if no other dislocations £j8 Xj , j ^ i were present in the crystal). As we will see, the self-energy 
of a single dislocation bS x is asymptotically (as e — > 0) independent of the position of the dislocation 
point x in f2, depending only on the elasticity tensor C and on the Burgers vector b. 

We begin by computing the self-energy in a very simple situation. We assume for the time being 
that Q is the ball B\ of radius one and center zero, that fi e = fx — b5o, with \b\ = 1, and we compute 
the energy of an admissible strain (3 considering the toy energy given by the square of the L 2 norm of 
0. More precisely, we consider 

(9) El°y{^P) := [ \/3\ 2 dx, 

JB 1 \B E 

and the self-energy induced by [i given by minimizing this energy among all admissible strains (3 
compatible with /j, s , i.e., satisfying 

(3-tds = b. 



In this setting the minimum strain /?(//) can be computed explicitly, and in polar coordinates it is 
given by the expression 

0{tx){6,r) :=-Lt(r,9)(g>b, 
2irr 

where f(r, 6*) is the tangent unit vector to dB r (0) at the point (r, ff). We deduce the following expression 
for the self-energy 



(10) 



£f f == / ^- dr = i(log 1 - loge). 
J s 2nr ztt 



Hence, as e — > the self-energy of a single dislocation behaves like | loge|. 

Notice that most of the self-energy is concentrated around a small region surrounding the disloca- 
tion. Indeed, fix s > 0, and compute the energy stored in the ball B E s(Q), 

f |/3( M )| 2 dx= f ^-dr= ^(log( £ *) - logs) = ^(1 - s)\\oge\ = (1 - s)Ef f . 

JB E s(0) Je 27Tr 2lT 27r 

As we will see, in view of Korn inequality, the logarithmic behavior and the concentration phenomenon 
of the self-energy hold true also in the case of clastic energies depending only on /? sym like in ([7|) . In 
this respect the position of the dislocation and the shape of £1 itself do not have a big impact on the 
value of the self-energy. It seems then convenient to introduce the self-energy as a quantity depending 
only on the Burgers vector b (and the elasticity tensor C) through a cell-problem. Before doing that, 
we proceed heuristically, by introducing the notion of hard core of the dislocation Sq, as a region 
surrounding the dislocation such that 

i) The hard core region contains almost all the self-energy; 

ii) The hard core region shrinks at the dislocation point as the atomic scale e tends to zero. 

To this aim, in view of the previous discussion, it is enough to define the core region as B Pe (0), where 
the hard core radius p e satisfies 

i) lim^o ff = oo for every fixed < s < 1; 

ii) p 2 — > as e — > . 

Note that condition i) is indeed equivalent to 
i') Um _ MM = 0. 

A direct consequence of our definitions is that the self-energy can be identified (up to lower order 
terms) with the elastic energy stored in the hard core region, for which, with a little abuse, we will 
use the same notation, i.e., 

(11) Ef { :=[ \(3(p)\ 2 dx. 
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Consider next the case of a generic configuration of dislocations p e := Xa=i b5 Xt in ^ (but for sake 
of simplicity we will keep the Burgers vector equal to b and the toy elastic energy (J9j) ) . The hard core 
region of \x e is given by the union of the balls B Ps (xi). 

Again we require a decay for p e and that the area of the hard core region tends to zero as e — > 0, 
N e — > oo, i.e., 

i) lim e ^o fr — o° for every fixed < s < 1; 

ii) N e p 2 £ -> as e -> 0. 

Moreover, we require that the dislocations are separated by a distance 2p e , i.e., 

iii) the balls B Pe (xi) are pairwise disjoint. 

This condition motivates the name of the hard core region. In view of this assumption, it is natural 
to identify the self-energy of p, e as the elastic energy stored in the hard core region. Therefore, we set 

Ef\p e ) := / |/% £ )| 2 dx = J2 l/3(Me)| 2 dx. 

J Hard Core i=1 J B Pe (a*) 

We expect each term of the sum in the right hand side to be asymptotically equivalent to the self- 
energy of a single dislocation introduced in (fTTj) . so that in view of (JTDJ) , 

(12) Ef V) = E / l^)! 2 dx K N *\ l0 S £ l- 

This expression represents the asymptotic behavior of the self-energy as e — > 0, N e — > oo. 

2.2. Long range interaction between dislocations. Since the self-energy is concentrated in a 
small region surrounding the dislocations, it is natural to define the interaction energy E E ntC1 ' as the 
energy diffused in the remaining part of the domain, far from each dislocation. In view of this definition 
we have that the total energy E e (p, e , /?(/z E )) is given by the sum of the self and the interaction energy: 

E e (He,P(Me))= [ W(f3(p £ ))dx= f W(/3(p e )dx 

jQ, £ (fi E ) ./Hard Core 

+ / W(0O*s)) dx = Ef l (p E ) + E™ ter ( Me ) , 

7n\Hard Core 

where we recall that f2 e (/z e ) is the region of outside the dislocation cores, i.e., O e (/i E ) := fi \ 
U±?iB s (xi), for fi £ = $xib- Let us compute heuristically the interaction energy for a very simple 

configuration of dislocations and for the toy energy considered in ([9|) . Let Q, be the unit ball B\ , and 
H e be a configuration of N £ periodically distributed dislocations whose Burgers vector b has modulus 
one. The stored interaction energy is given by 

/> 1*1 />27rr 

E e ntc '\^,f3(^ £ )):= \p(p £ )\ 2 dx = dr Xe |/3(r, 9)\ 2 d0 , 

JSl\Hard Core JO J 

where \e denotes the characteristic function of the set \ Hard Core. Thanks to the uniform distribu- 
tion of the Xi 's in B\ = Q we deduce that the number of dislocations contained in each ball of radius r 
is proportional to the area of such a ball, and more precisely is of order wr 2 N £ . Therefore, the average 
P t (r) of the tangential component of the strain on each circle dB r outside the dislocation hard cores, 
i.e., of f3t(r, 8)xe, is of order Trr 2 N e /2nr — N e r/2. The error is small thanks to the fact that the area 
of the hard cores region is negligible. By Jensen's inequality we obtain the following estimate from 
above for the interaction energy 

^ nter (^ e )= / / Xe\P{r,0)\ 2 d9dr> 2nr\P t {r)\ 2 dr w / 2nN 2 r 3 /Mr = CN 2 , 
Jo Jo Jo Jo 

where C is a constant independent of e. 

Note that the self-energy has been estimated looking at the circulation condition that the strain 

f3 e (fi) has to satisfy (in order to be an admissible strain) on the circles dB r (xi) for e < r < p e , while 

the estimate for the interaction energy has been established looking at the circulation condition on 

all circles dB r for r < 1. We will see that this estimate is indeed sharp for recovery sequences in 

T-convergence. We will then conclude that the interaction energy behaves like N 2 as e — > 0, N E — > oo. 
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2.3. Energy regimes. In view of the heuristic arguments of the previous section we fix a function 
s i— > N e that represents the number of dislocations present in the crystal corresponding to the internal 
scale e. The above considerations can be summarized in three regimes for the behavior of N s with 
respect to s — > 0: 

1. Dilute dislocations (N £ « |loge|): in this regime we have that the self-energy, which is of 
order N £ \ loggj, is predominant with respect to the interaction energy. 

2. Critical regime (N £ ss | logeQ: in this regime we have that the self-energy and the interaction 
energy are both of order N £ \ log e| « | loge| 2 . 

3. Super-critical regime (N £ y-> | logeQ: in this regime the interaction energy, which is of order 
|7V e | 2 , is predominant with respect to the self-energy. 

3. The setting of the problem 

In this section we specify the mathematical setting of the problem. In particular, we introduce the 
class X £ of admissible configurations of dislocations, the class AS £ (p) of the corresponding admissible 
strains, and the rescaled energy functionals T £ . 

From now on SI is a bounded open subset of R 2 with Lipschitz continuous boundary, representing 
a horizontal section of an infinite cylindrical crystal. For the given crystal, we introduce the class of 
Burgers vectors S — {b±, b s }. In what follows we assume that S contains at least two (independent) 
vectors, so that 

(13) Span R S* = R 2 

and this will imply that the function ip in the energy (JSj> is finite in whole of R 2 . The case of only one 
Burgers vector is easier and it implies that p is finite only on a one dimensional subspace of R 2 . We 
denote by § the span of S with integer coefficients (S = Span^S 1 ), i.e., the set of Burgers vectors for 
"multiple dislocations" . 

As in Section [2l N £ represents the number of dislocations present in the crystal, corresponding 
to the internal scale e. We introduce also the sequence p £ representing the radius of the hard core 
surrounding the dislocations, and we require 

i) lim e ^o fr = oo for every fixed < s < 1; 

ii) lim^ol^elp 2 = 0. 

Condition i) says that the hard core region contains almost all the self-energy, while condition ii) says 
that the area of the hard core region tends to zero, and hence that its complement contains almost 
all the interaction energy. Note that conditions i) and ii) are compatible whenever 



(14) N £ e s — > for every fixed s > 0. 

We assume that the distance between any pair of dislocation points is at least 2p £ and we define 
the class X £ of admissible dislocations by 

M 

(15) X £ := { M eVW(ft;R 2 ) : /x = ^ & S Xi , M e N, B Pe (a*) C SI, 

i=l 

\xj -Xk\> 2p e for every j ^ k, & G § j , 

where A4(Sl;M. 2 ) denotes the space of vector valued Radon measures. Given p, e X £ and r £ R, we 
define 

(16) Slr(ji):=Q\ (J B r (an), 

X;£ supp(/i) 

where B r (xi) denotes the open ball of center x\ and radius r. 

The class of admissible strains associated with any fi := X^i=i 6t ^ i s then 

\ j3 E L 2 {Sl £ (n):M 2x2 ) : Curl/3 = in Sl £ (p) and / f3-tds = ^ for all i = 1, ...,M \ . 

I JdB e { Xi ) I 
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Here t denotes the tangent to dB e (xi) and the integrand (3 • t is intended in the sense of traces (see 
Theorem 2, page 204 in [8] and Remark[T]for more details). In the arguments below it will be useful to 
extend the admissible strains to the whole of f2. There are various extensions that can be considered 
and that are compatible with the model that we have in mind. We decide to extend the /3's setting 
their value to be zero in the dislocation cores. Thus, from now on the class AS £ (p) of admissible 
strains associated with any jj, := ^Z i=1 & S Xi is given by 

(17) i5 £ M :={/?eI 2 (!l;M 2x2 ): [3 = in O \ O e (/i), Curl/3 = in O e (/*), 

/ (3-tds = £i, and / (/? - /3 T ) da; = for alH = 1, M }. 

JdB e {xi) J 

In view of the definition of the elastic energy, the last condition in the definition of AS e (p) is not 
restrictive and it is there to guarantee uniqueness of the minimizing strain. 

Remark 1. Let (3 E AS e (fi). Formally Curl/3 is a tensor defined by 

(0^1^,:=^.-^!, 

which by definition is antisymmetric with respect to the entries j, I. Therefore curl/3^) (where /3/j) 

denotes the i th row of /3) can be identified with the scalar 

d d 
curl/3 (l) := - — ft 2 - -5 — Pil- 
v ' 0x1 0x2 

In the sense of distribution, curl (3u\ is given by 



(18) <curl/3 ( i),^>= / < /3 (i) , JVcp >, 



where J is the clockwise rotation of 90°. From equation (|18p it turns out that whenever (3u-\ is in 
L 2 (fl;M. 2 ), then curl/3(j) is well defined for ip £ H 1 ^) and acts continuously on it; so that 

curl/% G H-^tyM 2 ) for every (3 € AS e (/x). 

On the other hand, if /3 (l) € ^(ftjIR 2 ) and curl/3 (l) e ff-^K 2 ) then /3 (i) is in L 2 (ft;IR 2 ) modulo 
gradients. 

Finally, notice that for every /i := X)f=i £* ^ -^ej the circulation condition in (|17p can be stated 
in the following equivalent way: 



M 



< Cwl /3,<p>= ^2 & 

»=i 

for every 93 6 #o(^) sucn * na * V = c * m Be( x i)- l n particular, if ip belongs also to Cg(fi), we have 

< Curl/3, </3 > = / <pd/x. 



S2 



The elastic energy E s corresponding to a pair (/i,/3), with /i£l E and /3 £ AS £ (/i), is defined by 



(19) E e (p,p) := / ^(/3)dx, 
where 

W(0 : : t 

is the strain energy density, and the elasticity tensor C satisfies 

(20) ci|f ym | 2 < C£ : £ < c 2 |f ym | 2 for any £ e M 2x2 , 

where £ sym := + £ T ) and ci and C2 are two given positive constants. Since /3 is always extended 
to zero outside J7 E (/x), we can rewrite the energy as 

E e (p,0)= f W (13) dir. 
Jn 
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Note that for any given \i G X e , the problem 

min / W(0) dx 

has a unique solution. The proof of this fact can be obtained by taking of a finite number of segments 
by fl e (/i) , obtaining in such a way a simply connected domain where to apply Korn inequality to the 
curl-free admissible strains, and then following the direct method of calculus of variations. 

We conclude giving the notation of the rescaled functionals in the different regimes for N S) according 
to the asymptotic analysis given in the previous section: 

Sub-critical or dilute regime: N e «. |loge|. In this case the predominant contribution comes 

from the self-energy and is of order N e \ loge|. Therefore, we set 

(21) T dn ^(u 0) - fw £ ^'^ if ^ G X *> 6 - 45e ^' 

k ' £ 1 \oo otherwise in L 2 (ft;M 2x2 ). 

Critical regime: N £ — |loge|. In this case both energies are of order A^|loge| — |loge| 2 . 
Therefore, we set 



(22) ?e{lhP) 




?E e {n,d) if/ieX e , 0e AS s {n), 

otherwise in L 2 (Q;M 2x2 ) 



Super- critical regime: N e » |loge|. In this case the interaction energy is the leading term, 
and is of order N 2 . Therefore, we set 

1 ' £ K \oo otherwise in L 2 (f7;M 2x2 ) 

Remark 2. The choice of representing the Curl constraint by defining the measure [i as a sum of 
Dirac masses concentrated in the dislocation points is one of the possible choices. Other possibilities 
would be to consider more regular measures with the same mass, such as 

M 

(24) fi = ^X*&Lt( Xi) 

%— 1 

or 

The advantage of these alternative choices is that jl and fi belong to iJ _1 (J7;R 2 ), which is the natural 
space for Curl/3. Indeed, since ft, (ft respectively) are in _ff -1 (£!;R 2 ), we could rewrite the class of 
admissible strains as follows: 

{0 e L 1 (f7;M 2><2 ) : Curl/3 = ft (respectively ft) in Q} . 

This notion of admissible strains does not coincide with that given in (|17p . but it turns out to be 
equivalent to (|17p . in terms of T-convergence, in the study of the asymptotic behavior of the energy 
functionals as e — > 0. 

4. Cell formula for the self-energy 

The self-energy stored in a neighborhood of a dislocation is of order |loge| and, in view of the 
concentration of the energy, it is asymptotically not affected by the shape of the domain, depending 
only on the elasticity tensor C and on the Burgers vector b. It seems then natural to introduce 
rigorously the notion of self-energy through a cell problem. 

In the following we will consider the self-energy of any multiple Burgers vector £ S S. For conve- 
nience we will introduce all the quantities we need for a generic vector £ £ M 2 . 

For every (el 2 and for every < n < r2 S M, let 

(26) AS ri ,r 2 (0 ■= <0eL 2 (B r2 \B ri ;U 2x2 ) :Curl/3 = 0,y 0-tds = A, 
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where B r denotes the ball of radius r and center 0. 

We first note that the admissibility conditions above on a strain (3, assure an a priori bound from 
below for its energy. This is made precise by the following remark. 

Remark 3. Given < r\ < r-i and ^gl 2 , for every admissible configuration (3 £ AS rijr2 (£) we have 



L 



\f3 s ^\ 2 dx>c\e, 

IB T2 \B T1 

where the constant c depends only on r\ and r-i- 

Indeed, by introducing a cut with a segment L, the set {B r2 \ B ri ) \ L becomes simply connected. 
Since Curl/3 = in (B r2 \ B ri ) \ L, there exists a function u £ H 1 (B r . 2 \ B ri ;M. 2 ) such that Vu = (3 
in {B T2 \ B T1 ) \ L. By the classical Korn's inequality applied to u we obtain 

|/3- A\ 2 dx < C [ \(3 sym \ 2 dx, 

'B r2 \B ri JB,, 2 \B ri 

for some skew symmetric matrices A. Moreover, by the fact that (3 £ AS rit r 2 (r]n), we conclude 

2 



r 2 1 

\P-A\ 2 dx > / — 

B r2 \B ri Jri 2lr P 

1 2tt P 



{13- A) -ids 



dp 



r2 1 ,., 2 , l€! 2 log: 



f 13 -ids dp = / 2 ^-\^\ 2 dp 
JdB„ in 2 ?rp 



(28) 



Set now C £ := Bi\ B £ , and let ip £ : R 2 — > K be the function defined through the following cell 
problem 

(27) M0~r^r— r min / W(/3) cfa; for every (el 2 . 

By (I2U1) . it is easy to see that problem (|2T|) has a solution. We will denote by /3 £ (£) the (unique) 
solution of the cell problem (|2"T|) whose average is a symmetric matrix. Note that /3 e (£) satisfies the 
boundary value problem 

DivC/? e (£) =0 in C s , 
Cpefc) -v =0 on dC s . 

Moreover, we will denote by /3r2 (£) : R 2 — > R 2 the planar strain defined in all R 2 corresponding to 
the dislocation centered in with Burgers vector £ . The strain /3r2 (£) is of the type 

(29) /3 E 2(£)(r,0) = ir s (0), 

where the function T depends on the elastic properties of the crystal, namely on the elasticity tensor 
C, and /3r2(£) is a solution of the equation (we refer the reader to [3] for a detailed treatment of the 
subject) 

' Curl & a (0 =tfo inM 2 , 
DivC/3 K 2(£) =0 inM 2 . 

In the next lemma we will see that in the cell formula lj2"T|) we can assign a suitable boundary 
condition without affecting the asymptotic behavior of the energy ip £ (£) . 

Lemma 4. Let £ £ R 2 be fixed and let (3 £ ASi t£ (£,). Assume that (3 satisfies 
(31) l/3(a0| < | -j- /or ei>en/ x £ C £ and for some K > 0. 

T/ien </iere exists a function [3 £ AS\ t£ {^) such that 

1) (3 coincides with (3 in a neighborhood of dC £ ; 

2) / W0) dx< I W(0 E (£)) dx + C\£\ 2 , where C > depends only on K. 
Jc E Jc, 



(30) 
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Proof. As in Remark [3] we make the annulus C £ a simply connected domain by introducing a cut with 
a segment L and we denote it by C £ := C £ \ L. Since both (3 and /? e (£) are curl free in C £ , we have 

/3 e (f)=Vti, (3 = Vv for some u, v G H^C^R 2 ). 

We want to construct f3 as the gradient of a convex combination of u and v. To this purpose, let 
ipi : (e, 1/3) — > M be a piecewise affine function defined by 

for r G (e, 2s), 

1 for r G (3e, 1/3), 
tpi is linear in (2e, 3e) 




and let ^2 : [1/3, 1) 

'l for r = 1/3, 

for r G [2/3,1), 

y>2 is linear in (1/3,2/3). 
Define 

Cl := {x G C* e : 2e < |x| < 3e}, C e 2 := {1 £ C E : 1/3 < |x| < 2/3}. 

Finally, set 

Q( X ) ■= / V (vi(NI)K^ - c x (u)) + (1 - Vl (\x\))(v(x) - Cl (v))) for \x\ < 1/3, 
\v(^ 2 (|x|)(n( a: )- C2 (u)) + (l-^ 2 (|a ; |))(t;( a; )-c 2 ( U ))) for \x\ > 1/3, 

where Cj(/) denotes the mean value of the function / on C\. It is easy to check that by construction 
[3 belongs to ASi and coincides with (3 in a neighborhood of dC £ . It remains to prove property 
2). 

By construction we have that (3 coincides with /3 e (£) in the region 

Cl := {x€C £ :3e< \x\ < 1/3}. 

Therefore, to conclude the proof of property 2) it is enough to compute the energy stored in C E \ C| . 
In view of (|20|) . of the fact that ||Vyi||oo < C/e, and of the following Poincare inequality 

\u - c(u)\ 2 dx < cr 2 \Wu\ 2 dx, 

(B 3r \B 2r )\L J (B 3r \B 2r .)\L 

where r > and c(u) stands for the average of u over the domain, it is easy to check that property 2) 
holds provided the following estimates are established 

(32) / \Vu\ 2 dx < C|£| 2 , / \Vv\ 2 dx < C|£| 2 , 

Jc e \ci Jc E \ci 

for some C independent of e. 

It remains to prove ([32|). By ([3T|l it follows that 

(33) / \Vv\ 2 dx<C\H\ 2 ( 

Concerning u, in view of l|28p and (|3H)) . we have that Vit can be written as 

(34) Vu = /MO + V/i, 
with /i given by 



,2e ! 




/. + 


I/ 3 * d 7 







Vh{x) = - VG(x,y)Cfo(t)'v(v)dy for x G C £ . 

Jasi(o)uas e (o) 

Here G(cc, y) is the Green function corresponding to the elasticity tensor C, satisfying the equation 



-Div^CVzG^,?/) = S y I in 



l 2 



for every fixed y. It is well-known (see [2]) that 
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(35) |VG(x)| < C/\x\ for every iel 2 , 

for some constant C > 0. From ([34]) and (l35l) it follows as in (l33l) that 



|Vu| 2 dx < 2 / (|/3 R 2(0| 2 + |V/f ) dx < C|£| 2 , 

C E \Cf JCe\CI 



which concludes the proof of (|32|) and therefore of the lemma. □ 

From Lemma HI together with (f2"9"|) , we deduce the following corollary 
Corollary 5. There exists a constant C > suc/i i/ia£ for every (el 2 , 

(36) < 7T^7 I W(fo(€))dx < A(0 ' ( Kr 



|loge|y Ce |loge| 
In particular, as e — > i/ie functions ip e converge pointwise to the function ip : R 2 — > K, defined by 

(37) ^(0 := lim ^(0 = lim n / W(PM0)dx=[ W(T^))d6. 

More precisely, we have 

C|C| 2 



(38) l^-^(OI< |Ioge| . 

By means of a simple change of variable e — > p we have that the self-energy is indeed concentrated 
in a p-neighborhood of the dislocation points whenever |logp| « | log e | . The precise statement is 
given in the next proposition. 

Proposition 6. For every e > let p £ > be such that \ogp e / loge — > as e — > 0. Let V> £ : R 2 — > R 
6e defined through the following minimum problem 

(39) := — j-, min { f W((3) dx : G AS e , Pe (£)}. 



Then ip e = ipeO- + o(s)), where o(s) —> as £ — > uniformly with respect to £. In particular, ip e 
converge pointwise as e — > to £/ie function ip : R 2 — ► R given in Corollary 
Moreover, let G ASe,^ (£), £ € R 2 , 6e sucft iftai 

(40) |/3(x)| <if|£ 

/or some if G R, and let ip £ : R 2 — * R oe defined through the following minimum problem 

(41) 4(0 := — !-r min { f W{0) dx : G AS e , Pe (£), ■ t = ■ t on dB £ U dB p \. 

lOE £ I. ./ c \ o J 



TTien ^ e = t/> e (l + o(e)), where o(e) — > as £ — ► uniformly with respect to £. In particular, ip e 
converge pointwise as £ — > to £/ie function ip : R 2 — > R given in Corollary 

Remark 7. The error o(e) appearing in the expression of ?/> and ^ in Proposition ((6|) can be estimated 
as follows o(e) « logp £ /log£. 

We are now in a position to define the density <p : R 2 — » [0, oo) of the self-energy through the 
following relaxation procedure 

{N N ~j 

X>fc^(&) = = JV€N, A fc >0, &€SI. 

fc=i fc=i J 

It follows from the definition that the function (/? is positively 1-homogeneous and convex. Moreover, 
since ip(£) > C|£| 2 f° r some C > 0, (that can be checked by its very definition), the inf in (|42[) is 
indeed a minimum. 
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Remark 8. Note that if for every z\, . . . , z s G Z we have 

(43) V>(X^iM >X^(&i), 

\i=l / i=l 

then in the relaxation procedure given in (|42[) we can replace 5 with S. More precisely, the formula 
for ip reduces to 

(44) <p(£) := min j £ |A^(6,) : £ AA = £, 6j G 5 1 . 

Actually, condition (|4"3")) can be viewed as a condition that the class of Burgers vectors bi of the given 
crystal has to satisfy in order to contain all the dislocation's defects observed in the crystal. In other 
words, if a dislocation corresponding to a vector b := J2i=i z ^°i stores an energy smaller than that 
obtained by separating all its components bi, then b itself has to be considered as a Burgers vector 
of the crystal. A rigorous mathematical definition of the class of Burgers vectors corresponding to a 
given crystal could be to consider the set {bi, . . . b s } C § of all vectors satisfying ip(bi) = f(bi), where 
§ is the set of slips under which the crystal is invariant. The Burgers vectors defined in such a way 
would always satisfy property (|4"3]) . 



5. A KORN TYPE INEQUALITY FOR FIELDS WITH PRESCRIBED CURL 

Let u G J^ 1,2 (f2; R 2 ), (3 its gradient and /3 sym and /3 skew its decomposition in symmetric and anti- 
symmetric part. The classical Korn inequality asserts that if f3 skcw has zero mean value, then its L 2 
norm is controlled by the L 2 norm of /3 sym . We will show that in dimension two the same result holds 
true also for fields (3 that are not curl free, modulo an error depending actually on the mass of Curl (3. 
The result is a consequence of some estimates for elliptic systems with L x -data recently proved by 
Bourgain, Brezis and Van Schaftingen in [3], [5]. The precise statement is the following 

Theorem 9 (A Generalized Korn type inequality). There exists a constant C depending only on Q 
such that for every (3 G L 1 (fl;M 2x2 ) with 

Cmlf3 = f ieM(n-,R 2 ), [ ^ skcw = 0, 

Jn 

we have 

(45) f \(3 s ^\ 2 dx<c( f \f3 s n 2 dx + (\^m)) 2 ) . 
Proof. The condition Curl (3 = (/ii, /Z2) can be written in the following form 

/ Wlf w ) xi = hi + m, 

I (/?12™% 2 = h 2 + t i 2 , 

where hi G i/ _1 (il) are linear combinations of derivatives of entries of f3 sym . 

Since the field ((Pf^™)^! > (Pi2 CW )x 2 ) is cur l lree > we deduce that curl(^i,^ 2 ) = curl(— hi, —h%), or 
equivalently 

(46) div (-/x 2 ,/xi) = div (h 2 ,-hi). 

By O Lemma 3.3 and Remark 3.3] (see also [25] and [3]) we have that if / G L 1 (f2;K 2 ) is a vector 
field satisfying div / G H~ 2 ( fl), then / also belong to H~ 1 (n) and the following estimate also holds 

||/IU-i(n) < c(||div /||jj-2 ( n) + ||/||i,i(o)) ■ 

This result clearly extents by density to measures with bounded variations. Thus, by (|46|) . we have 
that div (— H2,ni) G H~ 2 (fl), and so we deduce that fj, belongs to i7 _1 (J7) and 

(47) IHIff-vn) < c(||div (/ia, -hi)\\ H - Ha) + \fj\(Sl)) < c(||/3 sym || L2(f , ;R2) + | M |(0)) . 
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Now let u £ Hq(£1;M. 2 ) be the solution of — Au = (— /U2)Mi) m ^ an d let £ be the 2x2 matrix 
defined by £ := JVu (i.e., the i* 71 row of £ is given by £j = (— (ui) X2 , (ui) Xl ), for i — 1, 2). By definition 
we have that Curl£ = \i. By ([47]) we then obtain 

(48) / |£| 2 = ||Vu||i, < c\\n\\ H - HQ) < ( [ \(3 S ^\ 2 dx + (M(fi)) : 

Jn \Jn 

Since the average of the anti-symmetric part of £ can be easily estimated with the L 2 norm of £, 
we can assume that (j48|) holds for a matrix £ such that £ skow has zero mean value and Curl£ = fi. 
Therefore, by the classical Korn inequality applied to [3 — £, which by construction is curl free, and in 
view of (1481) we conclude 



f \l3 skcw \ 2 dx<c( f \/3 skcw -C kcw \ 2 dx+ f |f kow | 2 da; 
Jn \Jn Jn , 

< c (J |/3 sym - f ym | 2 dx + J |f kcw | 2 dx^j 

\(3 sym \ 2 dx + f \i\ 2 dx\<c([ |/3 sym | 2 rfx+ 01(0))' 
Jn J \Jn 



□ 



6. The critical regime (N e « | log e | ) 



In this section we will study the asymptotic behavior of the rescaled energy functionals as the 
internal scale e — > 0, in the critical energy regime, namely with N £ — | log e\ . In terms of T-convergence, 
we consider the rescaled energy functionals T £ : M(fl:R 2 ) x L 1 (fi;M 2x2 ) —> M defined in ([22]) . 

According to the heuristic arguments above, in this regime we aspect the coexistence of the two 
effects, the interaction energy and the self-energy, so that the candidate for the T- limit T : M.(VL\ K 2 ) x 
L 2 (f2;M 2x2 ) — > K is defined by 



(49) H,, S )Jl WWdx + SA^ dW 

y oc 



if /i e H-i^M 2 ), Curl (3 = n; 
otherwise in L 2 (fl;M 2x2 ). 

Theorem 10. The following I '-convergence result holds. 

i) Compactness. Let e n — > and let {(/i n ,/3 n )} be a sequence in M(fl;M. 2 ) x L 2 (Q;M. 2x2 ) 
such that T £rl {fi ni [3 n ) < E for some positive constant E independent of n. Then there exist a 
subsequence of e n (not relabeled), a measure fi £ Jf _1 (f2; K 2 ), and a strain (3 £ L 2 (fl; M 2x2 ) 7 
with Curl/3 = \x, such that 

(50) 7-j— — rMn ~^ M mM(fl;R 2 ), 

(51) TT^—^n^P mL 2 (n;M 2x2 ). 

|log£„| 

ii) T-convergence. The functionals T £ Y -converge to T as e — > with respect to the convergence 
in (|50p . (|51[) . i.e., the following inequalities hold. 

T-liminf inequality: for every (n,(3) £ (M(fl; M 2 ) P\ i? _1 (f2; M. 2 )) x L 2 (n;M 2x2 ), with 
Curl/3 = [i, and for every sequence (fi £ , (3 £ ) £ X £ x L 2 (Q; M 2x2 ) satisfying ([5u]) and (ISTj) , 
we Ziawe 

liminf T £ {^i £ ,(3 £ ) > T{n,(3). 

s— *0 

T-limsup inequality: giwen (/i,/3) S (.M(f2; M 2 ) n i7 _1 (f2; M 2 )) x L 2 (ft;M 2x2 ), with 
Curl/3 = jti, t/iere exists (/i £ ,/3 £ ) € 7W(fi;R 2 ) x L 2 (ft;M 2x2 ) satisfying ^ and ([51) . 

limsup.F e (/x £ ,/? e ) < Tdi, (3). 

6^0 
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Remark 11. Consider the functionals £ e {pi) := min^ T E {p, (3), representing the elastic energy induced 
by the dislocation measure /i. By the T-convergence result stated in Theorcm[in]wc immediately deduce 
that the functionals £ £ T-converge (as e — > 0) to the functional £ : M(fl; M 2 )ni/ _1 (il; M 2 ) — > K defined 

by 

(52) £(jx) := min{^(//,^) : Curl/3 = fj,}. 

Therefore, the energy £ (p) induced by a distribution of dislocations \x in the critical regime is given 
by the sum of both an elastic and a plastic term. In particular, any distribution of dislocations in this 
regime induces a residual elastic distortion (i.e., if \i is not zero, then so is the corresponding strain (3 
that minimizes (|52]l and hence its elastic energy). 

6.1. Compactness. Let {(/i n , /3 n )} be a sequence in M(£l; R 2 )xL 2 (17; M 2x2 ) such that T Zn (/^„, j3 n ) < 
E for some positive constant E independent of n. We give the proof of the compactness property stated 
in Theorem [TU] in three steps. 

Step 1. Weak compactness of the reseated dislocation measures. 

We first show that the sequence {(1/| loge n |)/Lt n } is uniformly bounded in mass. Let /i„ = X)i=i £ 
with £i >n £ §; we claim that 

1 1 M " 

(53) T ,1^1(0)=,! ,5>i )n |<C. 

|log£„| |loge„|^ 

Indeed, 

M„ 

E > T En (/i„, (3 n ) = M - ., / W((3 n ) dx > V n - , 2 / W(J3 n ) dx , 



loge n | 2 Jn e(M „) ^ |log£«| 2 is, 



, (*«,») 

where we recall that /3„ = in f2 \ f2 e (//„). After a change of variables we deduce 



^ E n 1 |2 / W(^n(^,n + y)) dy. 
trf I log e™ | 2 Jb d „ 



Note that the functions ?/ — > (3 n { x i,n + y) belong to the class AS 6n . Pcn (Ci.n) defined in (|26|) . Therefore 
we have 



m„ 

(54) £>En ^/ W((3 n ( Xi +y))dy 



1 Mn , Mn ... 

> E (e,«) = E &»l 2 ^ ( In 

|loge„|^ l lo S £ "l~t Vl&>l 



where the function ^ En is defined in ([39]) . Let ^ be the function given by Corollary and let 2c := 
inf|j| =1 By Proposition [6] we deduce that for n large enough 4> €n {0 > c f°r every £ with |£| = 1. 
By <J54]) we obtain 

where the last inequality follows from the fact that £i, n G § = Span^S 1 , S 1 is a finite set, and hence 
|£i, n | are bounded away from zero. We conclude that (|53|) holds. 

.fjtep 2. Weak compactness of the reseated strains. 
In view of the coercivity condition ([3]) we have 

(55) C|loge„| 2 > C\\oge n \ 2 T £ (fi n ,l3 n ) > C [ W(/3 n )dx > f \f3 s ^ m \ 2 dx. 

The idea of the proof is to apply the generalized Korn inequality provided by Theorem [9] to (3 n to 
control the anti-symmetric part of (3 n . Note that the Curl of j3 n is clearly related to the dislocation 
measure fi n , whose mass is bounded by C| loge„| by Step 1. On the other hand, it is not clear that 



16 A. GARRONI, G. LEONI, AND M. PONSIGLIONE 

|Curl/3„| < C|loge„|. Therefore we proceed as follows. For every in the support set of fi n , set 
Ci : n '■= B2s n {xi,n) \ B £n (xi yn ) and consider the function Ki^ n : Ci. n — > M 2x2 defined by 

In \x — Xi\ 

where J is the clockwise rotation of 90°. It easy to show that the following estimate holds true 
(56) f \K hn \ 2 dx<c[ \(%T\ 3 dx, 

Indeed, it is straightforward to check that 

\K hn \ 2 dx = C\^ n \ 2 
and by a scaling argument, from Remark [3] we get 

\Pr\ 2 dx>c\^, n \ 2 



and hence (j56j) . 

By construction Curl(/3„ — K i n ) — in CV„, and so /?„ — K i n = Vv i n in Cj.„ for some Vi >n S 
^(d^R 2 ). Thus by ® we have 



|V< y n m | 2 dx<C (|$Tf + |*Q,„| 2 ) dx < C / |/3^ m | 2 dz , 

and hence, applying the standard Korn's inequality to u, in , we deduce that 

/ |V^, n - A i>n \ 2 dx<C [ \VvZT\ 2 dx<C [ |/3^ m | 2 dx, 

where Ai_ n is a suitable anti-symmetric matrix. By standard extension arguments there exists a func- 
tion Ui tn £ H 1 (B2e„{xi i n);M 2 ), such that Vui, n = Vt>i,„ — Ai, n in C^ n , and such that 



(57) 



/ \Vu hn \ 2 dx< f \Vv i>n - A it7l \ 2 dx < C f |CT 2 dx. 



Consider the field n : ft -> M 2x2 defined by 

Pn(x) if a; G fi 



In view of (15711 it follows 



En ! 



Vu i<n {x) + A ii7l if x e B en (x^ n ). 



I |$Tf = / |/3**f dz + W |V< y ri 2 < C / |/?rT dz < C| log £ „| 2 . 

JO 4 JB en (x it „) JO 

By construction we have |Curl/3„|(f2) = |/i„|(f2); therefore, we can apply Theorem [5] to /3 n , obtaining 



(58) 



\k -A n \ 2 dx<C Qf |/3rf <*c + (lM«|(^)) 2 ) < C| loge„| 2 , 



where A n is the average of the anti-symmetric part of /?„. Since the average of (3 n is a symmetric 
matrix, we have 

(59) f \p n \ 2 dx= [ \/3 n - A n \ 2 - \A n \ 2 dx < [ \p n -A n \ 2 dx. 

Finally, by and (p)9")) we conclude 

|/?„| 2 dx< / \Pn-A n \ 2 dx< f n -A n \ 2 dx < C|loge„| 2 , 

which gives the desired compactness property for /3„/| loge„| in L 2 (Cl; M 2x2 ). 
Step 3. /i belongs to id^^il; R 2 ) and Curl/3 = /i. 
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Let <p £ Cq(Q). It is easy to construct a sequence {<p n } C Hq(FI) converging to <p uniformly and 
strongly in Hq(Q) and satisfying the property 

ip n = ip(xi tn ) in B £n (xi <n ) for every in the support set of p, n . 

By Remark [1] we have 

If 1 

ipdp = hm f / p n dp n = lim r-: 1 < Cwlp n ,(pn > 

(1 n^ + oo I loge„| J n rw + oo|log£„| 

= lim — r / (3 n JV(p n dx — [ [3J\7<pdx = < Curl/3, tp >, 

»^+oo |loge„| J n ' J n 

from which we deduce the admissibility condition Curl/3 = p. Moreover, since by the previous step 
we have /3 € L 2 (f7;M 2x2 ), we deduce that p belongs to if" 1 (H ; R 2 ). 

6.2. T-liminf inequality. Here we prove the T-liminf inequality of Theorem [TDJ Let 

{pi, (3) € M(fl;R 2 ) x i 2 (fi;M 2x2 ) with Curl/3 = p, 
and let (p £ , (3 £ ) satisfy ([50)) and (JHH). In order to prove the T-liminf inequality 

(60) limmi T e {p £ ,l3 £ )>T(p,f3), 

it is enough to show that inequality holds for the self and the interaction energy separately. More 
precisely, we write the energy corresponding to (p £ ,(3 £ ) in the following way 

Fe{Ve,Pe) = / XuB Pe (x» c ) W(A>) dx + / Xn\uB Pe (x» c ) W(/3 e ) dx 
Jn Jn 

and we estimate the two terms separately. Set r\ £ :— „• By Proposition [6] we have 

(61) nki L w * >- L * (£) * s /„ * (&)' ,+ *» * • 

where o(e) — > as e — * 0. Since Span R S = R 2 , the convex 1-homogeneous function ip defined in (|4"2"]) 
is finite in R 2 , and so continuous. Thus, in view of Reshetnyak Theorem (see [22], Theorem 1.2]), we 
deduce 

(62) ^ /„ * (£) * W (To^ I * 

From (|6ip and (|62p we have that the T-liminf inequality holds for the self-energy, i.e., 

(63) Ihninf ^ ^ XuB ,. ( „„^( A ) * > ^ „ ( ) *|. 
Concerning the interaction energy, by semicontinuity we immediately deduce 

(64) liminf * £ xo\us Pe (x 4 , e ) W(/3 e ) rfs 



In the last inequality we used the fact that the number of atoms of p £ is bounded by | loge| and 
hence, since by assumption | \oge\p £ — ► 0, the function Xn\uB p (x, e ) converges strongly to 1. 
Inequality (|63p . together with inequality (|64|) . gives the T-liminf inequality. 
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6.3. T-limsup inequality. Here we prove the r-limsup inequality of Theorem [TOl We begin with a 
lemma that will be useful in the construction of recovery sequences also for different energy regimes 
under consideration. Given fj, e := Yli=x ^i,e5 Xi e £ X e and r e — > 0, we introduce the corresponding 
measures, diffused on balls of radius r £ and on circles of radius r e , respectively, defined by 

i=i £ i=i e 

For every x e .i in the support set of fj, e we define the functions K^":\ K^f : B Tr {x £ ^) — » M 2x2 as 
follows 

(66) kf/{x) := -ij^i ® J(i - x e ,t), K&(x) := ^ E ,i ® J *"^' , 

27T7-| ' 27T |x — 2r e ,i| 

where J is the clockwise rotation of 90°. Finally, we introduce the functions K£ e : SI — > R 2 , : 
SI -► R 2 , defined by 

M E M e 

(67) := Y, ^XBr. (* M ), ^ : = £ («..,) • 

i=l i=l 

Note that 

(68) Curl Kjf = pT e ' - A? , Curl = -f?' . 

Lemma 12. Let N e — > oo &e satisfying fl4]l. £ := ^fc=i Afc£fc £ §, A/c > 0, A := Xwt Afc, 

/i := ^ da; emd r £ := 1/ (2y/KN^). Then there exists a sequence of measures fj. e = 53fe=i CfcMe * n A" e , 

™£/i /ig o/ i/ie type <^ s i SMC ^ ^of B Te (x) C SI, |x — y\ > 2r e /or every x, y in the support set of 

fj, e , and such that 

(69) \Vs\/N e A A fe dx in M(Sl), 

(70) i n M(Q;R 2 ), inMSl;! 2 ), ^ - £ «n R 2 ), 

(71) jl" - ^' ~* ° ^onsZymtf-^R 2 ). 

Proof. First, we prove the lemma for M = 1 and \i = £ dx with £ e S. 

For /i := £ dx we cover R 2 with cubes of size 2r e , we plug a mass with weight £ in the center 
of all of such cubes which are contained in SI, and we set fj. e the measure obtained through this 
procedure. Let us prove (fTTj) . all the other properties following easily by the definition of fi £ . Since 
H — j^- converges weakly to zero in L 2 (fl; R 2 ), by the compact embedding of L 2 in H^ 1 we have that 
H - ^ -> in iJ-^SljR 2 ). The fact that A ^f'' -> strongly in ff-^fijR 2 ) follows directly by 

(|68p. since K^ E /N £ — > strongly in L 2 (S1; R 2 ) (that can be checked by a very simple estimate). 

The general case with £ S R 2 and M > 1 will follow by approximating /i = £dx with periodic 
locally constant measures with weight £& on sets with volume fraction A&/A. In each region where 
the approximating measure is constant we apply the construction above and then we take a diagonal 
sequence. □ 

We are in a position to prove the r-limsup inequality of Theorem 1101 We will proceed in several 
steps. 

Step 1. The case fi = £ dx. 

Given £ S R 2 and (3 £ L 2 (S!;M 2x2 ) with Curl/3 = £ dx, we will construct a recovery sequence 
C X e , (3 e £ AS e (/i e ), such that (/i e ,/3 e ) converges to (£ dx,(3) in the sense of ([50|) and (|5Tjl . 

(72) l imsup _L^ f W (p £ )dx< f (W(p)+<p(t))dx, 

e^O |k)g£| z J n ,/ n 



and satisfying the additional requirement that (/3 e /| log e| — j3) -t tends to zero strongly in H 2 (3S1). 
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Let ip be the self-energy density denned in l|42[). and let Xk > 0, ^ € § be such that £ = ^k£,k 
and 

M 

(73) ^(0 = £ A ^fe)- 

fc=i 

Consider the sequence fi e := Yli^i Cm^e* < given by Lemma [T2l with N e = | log e| . Note that, since 
A^p 2 — > 0, we have that r e ^>> p e , and so fi E £ X e . 

Since the function K^ s f defined in ([66]) belongs to AS E , Pe (£ e ,i) an d satisfies condition (|40|) , by 
Proposition [5] for every x £i i in the support set of [i e we can find a strain j3 e ^ : £1 — » M 2x2 such that 

l) /9 e ,ie^S e , p .(6,i), 

3) n^b / ^OM = + o(e)) where o(e) as e -> 0. 

Extend /3 £jJ ; to be K^f in B rE (x e ^) \ B Pe (x e ^) and zero in fl \ (B re (x e ^) \ B E (x e ^)), and set 



2) • t = l<\; ■ t on dB £ (x £<i ) U 9B Pe (x £>< ) 



>:;■ 



(74) 

i=l 

Then 

4) Curl/3 e = -A?+Ae: 
where and p, £ are defined in ([65]) . Finally, set 

(75) 4:=|loge|/3-^+4, 
where is defined according to ([6"?]) . By Lemma [T2] and ([6"5| 

Curl-^L^^) = ( M - + - -fi—) = ( M - -^-) - in H^^M 2 ). 

|loge| |loge| | log £ | |loge| |loge| 

Therefore, we can add to (3 £ a vanishing sequence . — > in L 2 (fi; M 2x2 ), obtaining the admissible 
strain 

(76) Pe--=(Pe + Re)Xn e (p c )- 

In order to prove that the pair (fi £ , (3 £ ) is the desired recovery sequence we have to check the following 
properties 

i) (3 £ converge to (3 in the sense of definition ([51]) ; 

ii) The pair ((3 £ ,[x £ ) is a recovery sequence, i.e., 

lim / W((3 £ )dx= [ (WW) + <p(£)) dx. 
e ^°Jn Jo, 

To prove i), we first note that since M e ~ | log e | , r e ~ l/y/e, we have that 



(/..) I log e | V E(AIe )\n Pe (^) I log e | 

which implies that (3 £ /\ loge| is concentrated on the hard core region. Then by Lemma [H \l4\/N £ A 
Afc for every k, and by property 3) we have 

(77) lim 1 / W%3 e )dig= lim 1 / W(/3 e )dx 
e-o |loge| 2 7 n e -o |loge| 2 J n \n Pe (^) 

e ^ I & ' fe=i fc=i 

In particular, we deduce that /3 e /| loge| is bounded in L 2 (f2;M 2x2 ). Since the L 2 norm of /3 e /| loge| 
is concentrating on the hard core region, we conclude that /3 e /|loge| converges weakly to zero in 
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L 2 (tt; M 2x2 ). On the other hand, one can check directly that j\\oge\ converges strongly to zero 
in L 2 (Q). Recalling that also R £ /\ loge| — > 0, by ([75| and (|76|) we conclude that i) holds. 
Next we prove ii), namely, that the pair (/3 e ,/n e ) is optimal in energy. We have 



lim 



1 



£^o | logep 
Since (3 £ /\\oge\ in J L 2 (S1; E 
1 



W{(3 £ )dx = lim 



1 



W{\\oge\(3 + p £ )dx. 



W{P e ) dx = lim 



n t -6|loge| 2 Jo 
2x2 ), taking into account also (|77|) . we conclude 
1 



e^O I log el 2 



W{\\oge\f3)dx+ / W{l3 e )dx \ = / (W(/3) + <p({)) dx 



n 



l^o | loge| 2 

Finally, by the Lipschitz continuity of <9f2, from ((75|) and (|76|1 we also deduce that (f3 e /\ \ogs\ — (3) ■ t 
tends to zero strongly in H~i(dVl). 

Step 2. The case /i :— X);=i XAi£l dx. 

In this step we proof the T-limsup inequality in the case of [i locally constant, i.e., of the type 

L 

(78) n:= ^2xA,£,idx, 

i=i 

where Ai are open subsets of O with Lipschitz continuous boundary and £j <E M 2x2 . The construction 
of the recovery sequence is based on classical localization arguments in T-convergence and takes 
advantage of the previous step. 

Let us set (3i :— /3\—Ai, and let fj,i lE , f3i, £ be the recovery sequence given by Step 1 applied with Q 
replaced by Ai, with /3 = /?;, /j, = dx. Finally let us set f3 £ : fl — > M 2x2 and p £ e M(fl; M 2 ) as follows 



if x € Ai 



Me 



E^ 1 



By construction we have /i £ € X £ . Moreover since 
1 



|loge| 



Curl/3 e Lfi e (^ e )|| H -i(n ; E2) 



I log e I 



■0 



fl~3(9A,) 



from the previous step it easily follows that 

1 



I logs) 



Curl&LQ £ (// e )^0 



strongly in _ff ^f^R 2 ). Therefore we can easily modify this sequence f3 £ by adding a vanishing per- 
turbation in order to obtain the desired recovery sequence p £ . 

Step 3. The general case. 

In this step we show how to construct a recovery sequence in the general case (namely for a general 
dislocation measure /j, € A4(£l;M. 2 )). 

Let (/i, 0) be given in the domain of the T- limit T . In view of the previous step and by standard 
density arguments in T-convergence, it is enough to construct sequences (/j, n ,Pn) with Curl/3 n = ji n 
and with //„ locally constant as in (178p such that 



(79) n -»/3 in L 2 (0;M 2x2 ), /x„ /i in M 2X2 ) and \/J, n \(n) -> \n\(Sl). 

Indeed, under these convergence assumptions we get the convergence of the corresponding energies, 



i.e., 
(80) 



lim T{0 n ,^r, 

71 — *- + OC 



By standard reflection arguments we can extend the strain f3 to a function (3 a defined in a neighborhood 
A of f2, such that Curl/3,4 = [ia is a measure on A and moreover |/^|(<9£1) = 0. 
Let ph be a sequence of mollifiers, and define 



fh ■= Pa* Ph\-^t, Qh ■= MA *p fc Lf2. 
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For h large enough these objects are well defined in CI and Curl/^ = 137,. Clearly 

(81) / h ^/3inL 2 (0;M 2x2 ), ghdx fi in M(Cl;R 2 ). 
Moreover, since |/i,i|(9f2) = 0, we have 

(82) \g h dx\{Sl)^\ii\(n). 

Next, we approximate every gh by locally constant functions, precisely, we consider a locally constant 
function gh t k such that 

(83) \g h ,k - 0h|£oo(fi;m=>) -> as k — ► cx>, and gh,k~ghdx = 0. 

Let Th,k be the solution of the following problem 

{Curl 77^ = gh,k - gh in CI, 
Div r Kk = in fi, 

rh,k ■ t = in 50. 

By standard elliptic estimates we have 

(85) |?"/t,fc|i,2(n;M 2 x 2 ) < C|<7h,fc - 5/i|L 2 (n ; R 2 )- 

Finally, we set fh.k '■= fh + ^h,k- By (f84|) we have Curl//^ = gr^fc- Moreover, by (|83]) . |85|) we have 

(86) / M -»/ h ini 2 (fi;M 2x2 ) as fc -> oo. 

By (|8Tj) . (|83| . ((86|) . using a diagonal argument we can find a sequence (/^„,/3„) satisfying ([79| . and 
therefore (|80|). 

7. The sub-critical case (iV e « | log e | ) 

In this section we study the asymptotic behavior of the energy functional E e defined in (fl9|) in 
the case of dilute dislocations, i.e., for N s « | log er| . In terms of T-convergence, it means that we 
rescale E e with a prefactor A^ £ |loge|, with N s « | log e| . As we discussed in Section [21 in this case 
the self-energy for minimizing sequences is predominant with respect to the interaction energy (see 
Remark [TJ} . 

In contrast to the critical case, we have that the prefactors of strains and dislocation measures in 
the sub-critical case are different. Indeed, the natural rescaling for the dislocation measures is given 
by N E . On the other hand, in order to catch the effect of the diffuse energy associated with a sequence 
(fi £ ,(3 £ ) with bounded energy we have to rescale the strains by (N E \\oge\) 1 ^ 2 . These two quantities 
clearly coincide only in the critical case N £ = | log e| . The effect of a different rescaling for strains 
and dislocation measures is that in the limit configuration /i and (3 are independent variables, i.e., 
the compatibility condition Curl/3 = /i disappears in this limit. Actually the admissible strains in the 
limit are always gradients, i.e., Curl/3 = 0. Heuristically this is a consequence of the fact that the total 
variation of Curl/3 e is of order N E , so that Curl (f3 e / (N s \ logej) 1 / 2 ) vanishes. 

The candidate T-limit of the functional JTf lute : M{Cl;R 2 ) x L 2 (Cl;M 2x2 ) defined in ([21]) is the 
functional ^" dllute defined by 

(87) ^ diluto ( M ,/3):= {J n - ./<> r V</|/< 




W[[3) dx+ I ip ( ^- I d\fi\ if Curl/3 = 0; 



otherwise in L (CI 



2VO. M2x2\ 



The precise T-convergence result is the following. 



Theorem 13. Let N s — > oo be such that N e /\loge\ — ► 0. Then the following T-convergence result 
holds. 

i) Compactness. Let e n — > and let {(/z„,/3 n )} be a sequence in M(Cl;R 2 ) x L 2 (Cl;M. 2x2 ) 
such that ^ lute (/x„, /3„) < E for some positive constant E independent of n. Then there exist 
H G M(Cl]R 2 ) and [3 € L 2 (Cl; M 2x2 ), with Curl/3 = 0, such that (up to a subsequence) 

(88) j^^^fi in M(Cl;R 2 ), 
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(89) in 

ii) T-convergence. XTie junctionals T e T '-converge as e — > 7 urat/i respect to the convergence in 
(I88|) and (|80[) . to £/ie functional J^ dllutc defined in Q87[) . More precisely, the following inequal- 
ities hold. 

T-liminf inequality: for every (/x, /3) £ .M(f2; M. 2 ) x L 2 (fi; M 2x2 ), wii/i Curl/? = ; and for 
every sequence (/u e , /3 e ) <E X e x L 2 (Q; M 2x2 ) satisfying (JSHJ) and (|8"9"]) , we /iawe 

liminf^ dilutc (^,/3 £ ) > ^ diluto (^,^); 

£—►0 

T-limsup inequality: given (/z, /3) e A4(Q; M 2 ) x L 2 (f2; M 2x2 ), with Curl/3 = 0, there exists 
ificfie) e X £ x L 2 (n;M 2x2 ) satisfying JSSJ and §§§ such that 

limsup^ dilutc ( Me ,/3 e ) < T dilutc (fi,(3). 

Remark 14. The independence of strains and dislocation measures in the r-limit is a consequence 
of the fact that, in the dilute regime, the interaction energy is a lower order term with respect to the 
self-energy. Indeed, by the T-convergence result stated in Theorem [13] we immediately deduce that 
the functionals f^iiutc^ d enne d by 

£f utc (^):= min ^ Uute (/i, (3), 

/3eAS E (n) 

T-converge (as e -> 0) to the functional £ dilutc : M(n;M. 2 ) -> R defined by 

The energy £ dlluto (/i) represents the energy stored in the crystal induced by the distribution of dislo- 
cations /i in the dilute regime, and it is given only by the self-energy. 

The proof of Theorem [13] follows the lines of the proof of Theorem [TO] For the reader convenience 
we sketch its main steps. 

Proof of Theorem [7ffl The compactness property of rescaled strains and dislocation measures stated 
in (|88| and ()89|) can be proved with minor changes as in the critical case. Let us prove that in this 
case Curl/3 = 0. Let ip S Cg(fi) and let {<p n } C Hq(H) be a sequence converging to cp uniformly and 
strongly in ffg(fi) and satisfying the property 

ip n = (p(xi iTl ) in B £n (xi, n ) for every Xi^ n in the support set of /x n . 

By Remark [1] we have 

< Curl/3, if >= lim — — — - < Curl/3„,<^ n > 

^-"JVe^llogenl 1 / 3 

r N ^ 1 ( a v N *'* ( a n 
= hm —pr — — / ip n a>„ = lim —pr / (pdfi = 

e n ^0 | logSnl 1 ^ N en J n E „^0 | log£„| 1 /2 J Q 

from which we deduce Curl [3 — 0. 

Concerning the T-convergence result, the proof of T-liminf inequality is identical to that of the 
critical case, so that we pass directly to the proof of the T-limsup inequality. 

As in the critical case, classical localization arguments reduce the problem to the case of /i constant. 
(Note that the density argument used in the critical case is even easier in the sub-critical case, since 
no admissibility condition Curl/3 = /i is required.) 

The proof of the T-limsup inequality reduces to find a sequence {/J. e } C X e , with l/N e /i e — ^ £dx in 
M{n- : R 2 ), and a sequence (3 e £ AS s (/x e ), with l/(JV e „| \oge n \) 1 / 2 [3 e (3 in L 2 (f2;M 2x2 ), such that 

(90) lim sup— -rl r f W{[3 £ )dx< f (W((3) + <p(t)) dx, 

e^O JV e „| 10g£„| J n J n 

and satisfying the additional requirement that (/3 e /| log e| — (3) ■ t tends to zero strongly in H~^(dCl). 
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Consider the sequence \i £ := 53i=i £e,iS Xc ( given by Lemma[12] Construct the functions f3 £ : ft — > M 2 
as in (El) . Then set 



(91) &:=(N e \loge\)*0-K£'+0 e 
where K£ e is defined according to (|67|) . By (f68|) 

Curl / 8 B L_fi e (/i e ) = -/i? 



By its definition the density of the measure jl r e B /(N e \ loge|) 2 tends to zero uniformly as e — > 0. 
Therefore, we can add to /3 e a sequence i? e in L 2 (Q; M 2x2 ), with i? e / (N £ \ loge|) 5 — * 0, obtaining the 
admissible strain 

(92) Pe:=(Pe + Re)xn e (^)- 

In order to prove that the pair (/i £ , (3 £ ) is the desired recovery sequence, we have to check the following 
properties 

i) f3 s converge to (3 in the sense of definition (|89"1) ; 

ii) The pair (f3 £ , p £ ) is a recovery sequence, i.e., 



lim / W(j3 £ )dx = / (W((3) + iptt)) dx. 

To prove i), recalling that by Lemma [T^l \fi^\/N £ Xk dx for every k, we have 

1 /■ 1 M M 

(93) lim — / W0 e ) dx = lim — ]T (nM&) = £ A fc V(&) = 

We deduce that (3 E /(N e \ logel) 1 / 2 is bounded in L 2 (f2;M 2x2 ). As in the critical case, the L 2 norm 
of /3 e /| loge| is concentrating on the hard core region, so that $ e /\ loge| converges weakly to zero in 
L 2 (f2;M 2x2 ). On the other hand one can check directly that Kg c /(N e \ loge|) 1/ ' 2 converges strongly to 
zero in L 2 (17), from which property i) follows. 

Concerning ii), we prove that the pair (/3 e ,/i e ) is optimal in energy. We have 



lim 777-n nr / W(p s )dx = ]im—— ^ / W(\ logel/3 + (3 £ ) dx. 

Since j3 s /\ loge| — 1 in L 2 (fl; M 2x2 ), taking into account also (|9"3"|) . we conclude 



(W(P)+^))dx. 

Finally, by the Lipschitz continuity of <9f2, from (|91| and (|92|) we also deduce that ((3 £ /{N £ \ logel) 1 / 2 — 
(3) ■ t tends to zero strongly in iJ~2(<9Sl). □ 

Remark 15. The case N £ < C has been considered in [Sj, where the asymptotic behavior of the elastic 
energy for a fixed distribution of dislocations /i e = /i is provided up to the second order, and in [21] , 
where the problem of T-convergence induced by screw dislocations is addressed (with E £ {(3) :— H/3H!), 
without any assumption involving the notion of hard core region (essentially with p £ « e). 

We could extend the result given by Theorem[13]to the case N £ < C, obtaining a T-limit which has 
still the form as in (l87|) . but with \i := J^fLi d ^x 4 > where ^ £ §, and with 9? : § — * E defined now by 

{N N ~\ 

X>(&) = = WeN, £ fe esl. 
A:=l fe=l J 
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8. The super-critical case (N £ » | log er|) 

In this section we study the asymptotic behavior of the energy functionals E £ defined in (fl9|) in 
the super-critical case, i.e., for N £ » | loge|. In terms of T-convergence, it means that we rescale E s 
by N 2 , obtaining the rescaled energy functionals ^"| u P er defined in (f^5|) . As we discussed in Section O 
in this case we have that the interaction energy for minimizing sequences is predominant with respect 
to the self-energy. 

The natural rescaling for the strains in this case is given by N B , but we don't have any control 
on the total variation of the dislocation measure. As a consequence we will get a limit energy ^" su P er 
defined on L 2 (il; M^) depending on strains given by 

(94) F supcl (f3 sym ) := [ W(f3 sym )dx if (3 sym G L 2 (fl; M 2 * 2 ), 

Jn 

where WI^ denotes the class of symmetric matrices M sym in M 2x2 . The precise T-convergence result 
is the following. 

Theorem 16. Let N £ be such that N e /\ loge| — > oo as e — > 0. Then the following I '-convergence result 
holds. 

i) Compactness. Let e n — > and let {(/x n , n )} be a sequence in R 2 ) x L 2 (0; M 2x2 ) such 
that J- £ uper (fj, n , f3 n ) < E for some positive constant E independent of n. Then there exists a 
strain /3 sym G L 2 (Q; M 2x2 ), such that (up to a subsequence) 

(95) -WT^ m ^ ^ Sym m L2 ^ ' 

ii) T-convergence. The functionals J^upcr T-converge as e — > 0, with respect to the convergence 
in (|95j) . to the functional J? su P cr defined in (|94j) . More precisely, the following inequalities 
hold. 

T-liminf inequality: for every /3 sym G L 2 (Sl;M 2x I ^) and for every sequence (fj, e ,f3 e ) G 
X £ x L 2 (Q;M 2x2 ) satisfying we have 

liminf^f lpcr (^,/3 £ ) > T supm '{(3 sym ); 

T-limsup inequality: given f3 synl G L 2 (f2;M 2x ^) i/iere exists (p £ ,(3 £ ) e X £ x £ 2 (f7;M 2x2 ) 
satisfying ()95|1 suc/i £/iai 

limsup^ supcr (^,/3 £ ) < .F supcr (/3 sym ). 

Proof of Theorem \16\ The compactness property is simply due to the usual apriori L 2 bound for the 
strains /3| ym , while the T-liminf inequality comes simply by lower semicontinuity. 

The main difference with respect to the previous energy regimes is in the proof of the T-limsup 
inequality. Again the strategy is to approximating a special class of limiting configurations and then 
to proceed by density, but in this case it will be more convenient to approximate the strains (3 with 
C 1 functions, so that their Curl's are continuous. 

Thus, fix j3 G L 2 (fl; M 2x2 ) such that Curl/? is a measure fj, of the type fj, — g{x) dx with g continuous 
and let us construct a sequence {/J. £ } C X e (fi) and a sequence f3 £ G AS £ (/i £ ), with (3 £ /N £ — j3 in 
L 2 (0; M 2x2 ), such that the following T-limsup inequality holds true 

(96) limsup-^ / W{[3 £ )dx< [ W{f3)dx. 

e^o N£ Jn Jn 

Arguing as in the proof of Temma 1121 it is easy to prove that there exist Cgl depending only on 
II.9||l«=(o ; r 2 ) an d a sequence of measures fi E := ^i,s^x ie € X £ , with |£ i)£ | < C, such that, setting 

r £ := C/ViVe, we have B re (xij) G ft, \xi, e — Xj yS \ > 2r £ for every Xi >e , Xj, e in the support set of /j, e , 
and, finally, 

^ ± n in M(fl; R 2 ), ^ - /x strongly in H' 1 ^ R 2 ), 
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where /t e is defined according to (|65|) . Consider the functions K£ E defined in (|67| and set f3 e := 
N E f3 + . By ([681) we deduce 

Curl^LQ £ (// e ) = (jm- ^1) - in ff-^R 2 ). 

Therefore, we can add to /3 e a vanishing sequence ^ — > in L 2 (f2; M 2x2 ), obtaining the admissible 
strain 

Ar ■= {Be +Re)Xn s {^)- 

In order to prove that the pair (/i e ,/3 e ) is the desired recovery sequence it is enough to observe that 
— > in L 2 (£l; M 2x2 ). Indeed, by construction we have M £ < CN e , and therefore 

□ 

Remark 17. Note that in the super-critical regime we can not have a compactness property for 
the antisymmetric part of the admissible strains [3 £ , and indeed it is easy to exhibit examples where 
||/3| ym /iV e || L2 ( n . M 2x2) < C and ||/?| kcw /iV e || i 2( r 2. M 2x2) — > oo. Note that, since we do not have any 
control on the mass of Curl/3 e , we can not apply Theorem [9] 
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